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ABSTRACT 


The velocity distribution of electrons along the symmetry axis is measured by measuring the farther end cap electrode current with 
negative bias voltage applied across thoriated tungsten filament and alternately through the acquisition of resonance absorption 
signal of trapped electrons with the storage voltage i.e., for axial motion in a quadrupole Penning trap are presented. Electrons 
emitted from a thoriated tungsten filament with a range of energies from 1.5eV to 6.5eV are trapped at pressure of 2 x 10°8 torr 
and magnetic field of 0.05T. In a quadrupole Penning trap thoriated tungsten filament is cathode, closer end cap acts as a grid and 
farther end cap is the anode. The farther end cap is biased with positive potential respect to the thoriated tungsten filament, the grid 
is held at zero potential. The first order derivative of the anode current with respect to the grid voltage is directly proportional to 
Electron energy distribution function of electrons. Alternately, the area under resonance absorption signal of trapped electrons with 
the storage voltage is plotted. The derivative of area under resonance absorption signal of trapped electrons with the storage 
voltage is also plotted. Then the derivative of area under resonance absorption signal of trapped electrons with the storage voltage 
being converted to energy from E = eV is plotted. Now the energy distribution function is transformed to velocity distribution 
function by converting energy into velocity using the relation v = /2E/m. Thus obtained plot is fitted to Gaussian distribution 
function which is one-dimensional version of Maxwell’s velocity distribution function and it is fitted very well. The Gaussian 
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distribution of velocity is obtained, with a maximum that coincides with the storage potential applied across the thoriated tungsten 


filament, which is same as that of the velocity distribution of the thermionically emitted electrons from the filament. 


Key words: velocity distribution function, quadrupole Penning trap, Maxwell's velocity distribution function, Gaussian distribution 
function 


1. INTRODUCTION 


Electrons are confined in Penning trap by the superposition of electric field and magnetic field. The electric field confines electrons 
in the axial direction through an electric potential minimum and the magnetic field applied along the symmetry axis of the trap 
confines the electrons in the radial direction. Therefore the trapped electrons have motion in axial direction due to electric field and 
radial motion due to magnetic field. The Maxwell-Boltzmann distribution law is applicable to identical, indistinguishable particles 
with any spin (Eg. molecules of gas, trapped electrons under local thermodynamic equilibrium) and the Fermi-Dirac distribution law 
is applicable to identical, indistinguishable particles with odd-half integral spin that obey Pauli’s exclusion principle (Eg. free 
electrons in a metal). If we assume the electrons inside the trap penetrating into thermal ensembles, then we find Maxwellian 
velocity distribution for electron-electron collision time ~ 10us under local thermodynamic equilibrium [1,2,3]. 

When the bias voltage applied across the thoriated tungsten filament is increased at a constant filament current, the energy of 
the electrons entering the Quadrupole Penning trap increases, then the flux of low-energy electrons that are trapped is reduced and 
trapping signal disappears gradually. However, if we increase the filament current at the constant bias voltage, the flux of low- 
energy electrons increases and the trapping signal gradually restored. The resonance absorption signals of detected electrons 
appear, as electrons reach farther end cap electrode of the trap and those electrons whose energies lower than the trap potential 
are confined. When the bias is increased gradually at the constant filament current there is a loss of signal and disappears. The 
signal is restored when filament current is increased at the constant bias voltage. The increase in number density with the increase in 
flux of trapped electrons causes the space charge induced effects and the signal vanishes completely, that indicates no trapping [4]. 
Even though the expected distribution function of thermionic electrons is Fermi-Dirac distribution [5], the distribution function 
obtained from the current measured on the farther end cap electrode results from the electrons that have travelled through the 
trapping region. The distribution function is expected to change due to the effects of imperfections in quadrupole geometry and 
collisions of electrons with the background gases in the trap. However, in these situations, thermalization processes, yielding a 
Maxwell distribution for electron-electron collision time ~ 10s, occur on a much shorter time scale [1] than it takes for the electrons 
to reach the farther end cap, for typical number densities of 101 m7? and energies of a few eV in the trap [6]. 

A graph of area under the detection signal is plotted versus the storage voltage. The area under the resonance absorption signal 
of electrons is directly proportional to the number of electrons trapped in ion trap. The derivative of area under the resonance 
absorption signal of electrons versus the storage voltage gives the energy distribution function. The data obtained is fitted to a 
Gaussian function and it is found to fit very well [3,6,7]. The energy is transformed to the velocity using the relation v = {2E/m. The 
graph of derivative of area under the resonance absorption signal of electrons versus velocity yields the velocity distribution of 
trapped electrons along the symmetry axis of the trap under axial oscillations. The velocity distribution is the Maxwell-Boltzmann 
distribution along the symmetry axis of the trap and in one-dimension which is a Gaussian distribution function. 

The velocity distribution function of trapped electrons obtained through the measurement of the energy of the axial oscillations 
in a quadrupole Penning trap. Electrons of energies ranging from 1.5eV to 6.5eV are emitted from a thoriated tungsten filament and 
are trapped at the pressure 2 x 10~® torr under magnetic field of 0.05T. Subsequently the trapped electrons are detected through 
monitoring their axial oscillations by an electronic tank circuit that is weakly coupled to the trap, by resonant energy transfer to the 
electrons. The Lab VIEW program technique developed enables us the direct measurement of the velocity distribution function of 
the electrons along the symmetry axis of the trap under axial oscillations. We obtain a normal distribution of velocity, with a 
maximum that coincides with the potential applied across the thoriated tungsten filament that thermionically emits electrons, 
indicates that the velocity distribution of electrons in the trap reflects the velocity distribution of the thermionically emitted electrons 
from the filament [3]. 


2. THEORY 


According to Maxwell, at a given temperature when the identical, indistinguishable particles are in thermal equilibrium, then the 
random motion can be grouped. The velocity distribution function f(v) gives the probability of velocities of particles that can 
occupy each state of velocity ranging from v and v + dv. When the initial minimum energy E;> kT, the Bose-Einstein and the Fermi- 
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Dirac distribution laws reduce to Maxwell-Boltzmann distribution law. The Maxwell's velocity distribution is skewed as it is not 


symmetric in velocity and given by [8]. 


_dny _( m 24 2 mv? : 
i ae = (ser) alle a Ne Ps ” 


For an electron the velocities at a temperature of 10000K are 


fox f 8kpT [ax iF 
= _ = 5.5 x 105ms~1,0 = =< = 6.2 x 105ms~!, Vins = = = 6.7 x 105ms~! (2) 
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The velocity v has components vx, vy and v, and the volume element is dv = dv,dvydv,. The kinetic energy of particles in all the 


three-directions is 
==mvz + =mvy +=mvz (3) 


When we introduce the functions $(v,), (vy) and o(v,) such that 
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The above particle velocity distributions are called Maxwell distributions. The average kinetic energy of particles along different 


directions is 


1 1 1 3 3 3 
5 mv?) = =m(v2) + =m(vz) + 5m(v2) = 5m(v2) = 5m(vz) = Smlv7) (6) 


Now, the average kinetic energy along x, y and z are equal according to the principal of equipartion of energy. 
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Thus, the particle kinetic energy per degree of freedom is ae and correspondingly the average kinetic energy of particle in the 
three-dimensional space is ye mv?) = ot Ts These relations define the temperature and hence average kinetic energy of particles 


[9]. The velocity distribution of electrons along the symmetry axis of the trap, which are under axial oscillations and velocity range 


7 m mv,? 3 
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The probability of finding an electron whose energy range between E, and E, + dE, along the symmetry axis of the trap is 


between v, and v, + dv, is a Gaussian distribution function given by [9] 


dP,(Ez) = ew [- (=) dE, (10) 


The probability density of an electron whose energy range between E, and E, + dE, along the symmetry axis of the trap is [10] 
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G)=— ( = ) 1 
pe(Es) = p—pexp (— Ea (11) 
The average energy of an electron along the symmetry axis of the trap is 
(E,) = i E,p,(E,)dE, = kpT © 1.38 x 107197 (12) 

0 

The form of three dimensional quadrupole electric potential is [11,12] 
= Yo 2 2 2 
Vi, y,zZ) = Pl —x*-y?*) (13) 


Fig. 1(i): Penning trap, B: Magnetic Field, Vo: storage voltage, (ii) Magnified view of motion of an electron in a Penning trap showing 
reduced cyclotron, magnetron and axial motions 


The vector potential of the magnetic field B is [9] 


1 1 Bot 4V, . 4V)B 
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The Lorenz force in an electromagnetic field with an electric field E and magnetic field B on electron of mass m, charge e, moving 


with a velocity v, is given by [9] 


F=el[E+(vxB)] (15) 
2Vo . 2V . 4VYo sia oes 
F= e |(e x8 + yy - 222) + BOR - 29) (16) 


The potential well due to the electric field in the axial direction causes harmonic oscillations of an electron with frequency given by 


1 |4eV, 
te 2m Ima? ca 


An electron in magnetic field performs a cyclotron oscillations in the absence of electric quadrupole field with cyclotron frequency 


[11,12] 


given by [9,11,12] 


l=5— (18) 
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The defocusing electric force modifies the oscillations of electron. The net force experienced by an electron in an electromagnetic 


field is balanced by the centrifugal force [9,11,12]. 


mv? _ 3 2eVo as 
- ev aD a (19) 
2 
w — ww += 0 (20) 
+ 2 2f2 
fy = Ba le fe i (21) 
The reduced cyclotron frequency is [9,11,12] 
Pal fe —2i7 
f= bOI the Z (22) 
This arises due to the oscillations of the electron in presence of bulk of electrons. The magnetron frequency is [9,11,12] 
—./f2 —2F2 
fn =e (23) 


2 


This corresponds to (£ x B) drift in the radial plane around the central region of the trap. 


Reduced Cyclotron 
/ Motion 


Resulting Motion 


Fig.2: Resulting motion of an electron due to the combination of reduced cyclotron, magnetron and axial motions 


Vo B re re fr re 
(V) (T) (MHz) (MHz) (MHz) (MHz) 
4V 0.05T 1399.164 1398.90675 0.5145 26.8288 
Table 1: Values of frequencies of an electron due to reduced cyclotron, magnetron and axial motions 
fini fri fe. © 1:52:2719 


3, EXPERIMENTAL PROCEDURE 


Measurement of velocity distribution of trapped electrons along the symmetry axis of the trap: The resistively heated 
thoriated tungsten filament emits electrons due to thermionic emission and all emitted electrons enter the trap volume through a 
small orifice of size 3mm in closer end cap, they can travel and reach the farther end-cap electrode of the trap [3,7]. The current on 
the farther end-cap is due to the electrons emitted by the filament irrespective of their kinetic energy. The integration of the 
electron energy distribution function f(E,) of thermionic electrons at farther end-cap represents all energy values. At lower bias 
voltages the electrons of low kinetic energy do not reach the end cap. The bias voltage is reduced till zero value; the current is 
obtained due to electrons with kinetic energy of an upper bound of the energy due to the bias voltage. The filament bias is varied 
between 1.5V to 6.5V that causes the entry of electron into the trap having different velocity distributions centered around the bias 
voltage are trapped. The energy of electrons is a few eV and the temperature of electrons is around 10000K which is much greater 
than the surrounding neutral gas molecules whose temperature is around 300K. The electron plasma generated is non-thermal 
plasma under local thermodynamic equilibrium. In this experiment electron-electron collision time tg, ~ 10us during which local 
thermal equilibrium is achieved [3], which is much smaller than the time (few tens of milliseconds) during which the measurements 
are carried out and hence we are probing the regime wherein the electron plasma has reached local thermodynamic equilibrium 
[3,9]. We find Maxwellian velocity distribution, for electron-electron collision time ~ 10us, when the electrons inside the trap 
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penetrating into thermal ensembles [1,3,9]. The graph of farther end cap electrode current versus negative bias voltage applied 


across the filament, and alternately the graph of area under detection signal of electrons versus storage potential is similar to the I-V 
curves obtained using Langmuir probes and yield energy distribution function and hence velocity distribution function. 


(i) Measurement of velocity distribution from farther end cap electrode current: In a quadrupole Penning trap thoriated 
tungsten filament is cathode, closer end cap acts as a grid and farther end cap is the anode. The farther end cap is biased with 
positive potential respect to the thoriated tungsten filament, the grid is held at zero potential. The farther end cap electrode current 
is measured with the variation of negative bias voltage from 1.5V to 6.5V which is applied across the thoriated tungsten filament to 
trap electrons. The first order derivative of the anode current with respect to the grid voltage is directly proportional to Electron 
energy distribution function of electrons [13,14]. The anode current J decreases with decrease in bias voltage from V to V — AV, then 
electrons whose energies E < eV do not reach the farther end cap electrode. Energy distribution function of electrons is obtained 
from the first order derivative of the anode current with respect to the grid voltage [3,15]. The reduction in current at the farther end 
cap is [15] 


co 


oN = él =e | f(Ez) dE, —e | f(E,) dE, (24) 
V (V-AV) 


aCe) [-w= [=e renee, (25) 
eldV edV # e J oe 


(ii) Measurement of velocity distribution from Resonance absorption signal strength: The detection of electrons is achieved by 
coupling the trap to a weakly excited detection circuit using a non-destructive detection technique. The confinement times of 
electrons in the quadrupole Penning trap in our experiments are in the range of few tens of milliseconds as shown in Fig.3, which is 
much greater than the time for the entry of the electrons into the trap to achieve the local thermal equilibrium which is around 10ys. 
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Fig.3: Signal fall-off with time with confinement time, t = 73.40 ms, at pressure P = 2 x 10-® torr, and magnetic field B = 0.05 T 


0.00 0.01 0.02 0.03 0.04 


The potential V is reduced to V' = V — AV, and then the electrons with energies E, < eV trapped only, then AN be the number of 
reduction in trapped electrons. The number of electrons trapped at the trap voltage Vis [3,9,15] 


e(V-AV) 


eV 
AN =e | f(E,) dE, -e i f(E,) aE, (26) 
0 0 
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Taylor expansion around eV, considering only first order term, neglecting higher order terms 


and taking limit as AV > 0 yields 


AN dN 


Aim av =a > f (Ez) (27) 


The derivative of the area under the resonance absorption signal of trapped electrons versus energy of electrons is plotted and 
slope of the curve gives the Maxwell-Boltzmann distribution function [3,9]. The graph yields the energy distribution and hence 
velocity distribution of trapped electrons along the symmetry axis of the trap as we measure only the axial motion [3,9]. The graph 
of the data is fitted to the following form of Gaussian distribution which is one-dimensional version of Maxwell-Boltzmann 
distribution for electrons trapped in ion trap [9]. 


A Xx — X_\?2 
y= yo + ern |-2 - )| (28) 


Electrons are loaded into the trap whose energy is distributed around the filament bias voltage, which is varied from 1.5 to 6.5 V 
that determines the energy and hence velocity of the electrons. Initially the storage potential is held at 6.5 V and then reduced to a 
potential V as shown in Fig. 4, which causes electrons with energy higher than V leaving the trap. After a certain time t, the detection 
of trapped electrons is carried out and the resulting signal strength is recorded. The measurements are made continuously by 
reducing the potential V in steps of 0.5 V, this voltage may be reduced to a minimum of 1.5 V. The loading of electrons into the trap 
volume is terminated by gradually reducing the filament bias to 1V. No detection is observed below this voltage as the storage 
potential cannot be lower than the detection potential, when no electrons are trapped. Then the electrons of all energies enter the 
trap, but those electrons whose energy is lower than the trap voltage are trapped. The potential is reduced to V, then the electrons 
whose energies greater than eV leave the trap. The time called the dwell time, during which the trap voltage remains V, then the 
electrons whose energies exceed eV escape from the trap. The dwell time is well within the confinement time of the trapped 
electrons. The voltage is then ramp down, and then the axial frequency of the trapped electrons is varied. At a certain ramp voltage 
the corresponding axial frequency matches the tuned detection circuit frequency, then the energy is transferred from the detection 
circuit to the electrons and the corresponding signal of electrons obtained is recorded. 


Step down 


Storaga voltage (V) 


T 7 T T T 
1.000 1.005 1.010 1.015 1.020 


time (Ss) 
Fig.4: Schematic diagram showing the variation in applied storage voltage with time. The electrons are injected in to the trap for 1s, 
and signal of electrons appears when the trap potential is ramped. 


When the detection circuit was tuned to an axial oscillation frequency that corresponded to about 1V, the resonant energy 
transfer takes place from the circuit to the trapped electrons; therefore the minimum storage voltage is set at 1.5V. Measurements 
on energy and hence velocity along symmetry axis of the trap are carried out for different storage voltages; the strength of 
resonance absorption signal of electrons is directly proportional to the number of electrons trapped and are recorded. 
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4. RESULTS AND ANALYSIS 


(i) Measurement of velocity distribution from farther end cap electrode current: The thermally generated electrons are 


accelerated, and then pass through the center of the trap and reach farther end cap electrode which results current. The negative 
bias voltage was varied from 1.5V to a maximum of about 6.5V in steps of 0.5V and the farther end cap electrode current was 
recorded at each bias voltage. All the electrons emitted by the process of thermionic emission from the thoriated tungsten filament 
can travel and reach the farther end cap electrode; the current is due to all the electrons being emitted irrespective of their kinetic 
energy. The farther end cap current is due to the energy distribution function f(E,) of the thermionic electrons over all energy 
values range between E, and E,+ dE, along the symmetry axis of the trap. At lower bias voltages the electrons which have low 
kinetic energy do not reach the farther end cap. The electrons at OV bias voltage, the kinetic energies have an upper bound of the 
energy due to the bias voltage. A graph of farther end cap electrode current versus the negative bias voltage applied across the 
thoriated tungsten filament is plotted, which shows that if there is no bias voltage on the farther end cap electrode there is zero 
current, when the bias is increased then the current increases gradually and reaches almost a constant value. A graph of farther end 
cap electrode current versus the negative bias voltage is plotted as shown in Fig. 5. 
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Fig.5: Farther end cap electrode current versus negative bias voltage for electrons 
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The farther end cap electrode current is converted into number of electrons using N = I/e and the negative bias voltage is 


converted into energy of electrons using E = eV. A graph of number of electrons versus energy of electrons is plotted, which shows 
the energy distribution of thermal electrons as shown in Fig. 6. The energy of electrons is converted into velocity using v = J2E/m 
and a graph of number of electrons versus velocity of electrons is plotted, which shows the velocity distribution of thermal electrons 
as shown in Fig. 7. 
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Fig.7: Velocity distribution of electrons measured from farther end cap electrode current 


(ii) Measurement of velocity distribution from Resonance absorption signal strength: The detection signal of electrons appears 
as soon as electrons enter the trap; the electrons which have energies lower than the trap potential are confined. When the bias is 
gradually increased, at the constant filament current, there is a loss of signal. At constant bias the signal is restored when the 
filament current is gradually increased. The restored signal is now at a different trap voltage on the ramp, when it is compared to the 
previous measurement. The trapping voltage shifts significantly with the variation of bias, due to the variation of the number of 
trapped electrons, this indicates the effect of space charge in shifting the trapping potential [16]. The appearances and 
disappearances of the signal alternately, when the bias voltage is varied is due to the velocity distribution of the electrons which are 
entering the trap. Thermalization processes that yield a Maxwell distribution due to electron-electron collisions, for typical number 
densities of 1015 m~3 and energies of a few eV in the trap which occurs on a much longer time scale than it takes for the electrons to 
reach the farther end cap. The axial oscillation frequency of the electrons is modified by space charge effects of cloud of electrons. 
The signal of trapped electrons is acquired and the area under the signal is found for different storage voltages. The area under the 
resonance absorption signal of trapped electrons is directly proportional to the total number of trapped electrons [3,9]. The number 
of electrons at the centre of the trap is measured using non-destructive detection technique at different trapping potentials. The 
axial oscillatory motion of the electrons is controlled through the Lab VIEW program that controls storage voltage. The number of 
electrons at the centre of the trap depends upon the storage potential. A graph of area under the detection signal is plotted against 
the storage voltage as shown in Fig. 8. The graph shows that the area under detection signal increases gradually and has plateau 
region beyond 6V. 

The derivative of the area under the resonance absorption signal of trapped electrons is plotted against the storage voltage and 
the function is fitted to a Gaussian function as shown in Fig.9 [3,9,17]. The storage voltage is converted into energy of electrons 
using E = eV and the graph of the derivative of the area under the resonance absorption signal of trapped electrons is plotted 
against the energy of electrons and the function is fitted to a Gaussian function which shows the energy distribution function f (E,) 
measured along the symmetry axis of the trap at the pressure of 2 x 1078 torr and the magnetic field of 0.05T and as shown in 
Fig.10. 
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Fig.8: Variation of area under signal with storage voltage at the pressure of 2 x 1078 torr and magnetic field of 0.05T 
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Fig.9: Variation of derivative of the area under the detection signals versus storage voltage at the pressure 2 x 10~® torr and 
magnetic field of 0.05T 


The energy of electrons is converted into velocity using v = J2E/m and a graph of derivative of the area under the resonance 
absorption signal of trapped electrons versus velocity of electrons is plotted, which shows the velocity distribution of thermal 
electrons f(v,) measured along the symmetry axis of the trap at the pressure of 2 x 10-8 torr and the magnetic field of 0.05T and as 
shown in Fig.11 [3,9]. 

Fig. 5 is the graph of area under the detection signal versus storage voltage and Fig. 6 is the corresponding derivative of the area 
under the detection signal versus storage voltage, which shows the energy distribution function. 
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Fig.10: Variation of derivative of the area under the detection signals versus energy at the pressure 2 x 10~® torr and magnetic field 
of 0.05T 
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Fig.11: Variation of derivative of the area under the detection signals versus velocity at the pressure 2 x 1078 torr and magnetic field 
of 0.05T 


5. CONCLUSION 


The trapped electrons that form non neutral, non thermal plasma are in local thermodynamic equilibrium under low pressure of 
~2 x 1078 torr and magnetic field of 0.05T. The strength of the detection signal of electrons varies with the storage voltage and this 
variation depends on the range of velocities of electrons due to the bias voltages applied across the filament. However, the variation 
of detection signal of electrons along the symmetry axis of the Quadrupole Penning trap is independent of the applied external 
magnetic field. This demonstrates the possibility of the direct measurement of the distribution function of the velocity along the 
symmetry axis of the trap and it is monitored through detection of resonance absorption signal in non-destructive electronic 
detection technique. The axial oscillatory motion of trapped electrons along the symmetry axis of the Quadrupole Penning trap is 
decoupled from their radial motion, under local thermodynamic equilibrium condition; therefore the velocity distribution function of 
the trapped electrons can be measured along the symmetry axis of the trap. It was found that the most probable energy of electrons 
is 4eV which is approximately 50% of the potential depth of the trap. The velocity distribution function shows that the maximum is 
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at the storage voltage that corresponds to the bias voltage across the thoriated tungsten filament. Thus, the distribution function of 


the electrons reflects the total velocity distribution measured along the symmetry axis of the trap, independent of the fact that 


whether we monitor only the axial motion, due to the consequence of equipartition of energy in all degrees of freedom of motion. 
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